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Context. The determination of the global structure of the planar and extra-planar Milky Way Hi disk depends critically on a reliable 
database but also on reasonable assumptions about the shape of the Milky Way rotation curve. 
Aims. We derive the 3-D Hi volume density distribution for the Galactic disk out to ~ 60 kpc. 
. Methods. Our analysis is based on parameters for the warp and rotation curve derived previously. The data are taken from the 

Leiden/ Argentine/Bonn all sky 21 -cm line survey. 
' (— I ' ' Results. The Milky Way H i disk is significantly warped but shows a coherent structure out to ~ 35 kpc. The radial surface density 

distribution, the densities in the middle of the warped plane, and the H i scale heights all follow exponential relations. The radial scale 
I . length for the surface density distribution of the Hi disk is 3.75 kpc. Gas at the outskirts for 40 :s ^ 60 kpc is described best by a 

O ' distribution with an exponential radial scale length of 7.5 kpc and a velocity dispersion of 74 km s"' . Such a highly turbulent medium 

^ ' fits also well with the average shape of the high velocity profile wings observed at high latitudes. The turbulent pressure gradient of 

C/3 , such extra-planar gas is on average in balance with the gravitational forces. About 10% of the Milky Way Hi gas is in this state. The 

. large scale Hi distribution is lopsided; for R 15 kpc there is more gas in the south. The Hi flaring indicates that this asymmetry is 

caused by a dark matter wake, located at ~ 25 kpc in direction of the Magellanic System. 

Conclusions. The H i disk is made up of two major components. Most prominent is the normal H i disk which can be traced to /? ~ 35 
kpc. This is surrounded by a patchy distribution of highly turbulent gas reaching large scale heights but also large radial distances. 
At the position of the Sun the exponential scale height in the z direction is 3.9 kpc. This component resembles the anomalous gas 
discovered previously in some galaxies. 

00 ' Key words. Galaxy: disk - Galaxy: halo - Galaxy: structure ~ Galaxy: kinematics and dynamics - ISM: structure 

' 1 . Introduction the gas within the disk we may proceed to determine the nature 

OO of the gas around this disk. 

The gaseous content of galaxies is a crucial counterpart to the 
^ stellar population but has the additional advantage that it can be 
• I— I traced far beyond the stellar population, radially but also per- 
pendicular to the disk. Stars are born out of gas condensations 
H and their deaths affect the three-dimensional evolution of the 

Galactic disc. The gas traces such activities. Here we are inter- In Sect. 2 we briefly introduce the database and the conver- 
ested in the steady st ate situa tion caused by dynamical processes sion of brightness temperatures Tb(1, b, v) to densities n(R, z, tp). 
as described e.g. bv Ide Av illez (2000). The 21-cm line traces The observer's location inside the Galactic disk causes some am- 
these processes, in particular the diffuse part of the interstellar biguities and constraints that need to be discussed. In Sect. 3 we 
medium (ISM). A thin disk is expected in the plane, overlaid derive radial dependencies of the average surface density dis- 
by a thick H i disk which is associated with an ionized layer, tribution within the disk. This is foUowed in Sect. 4 by a dis- 
Depending on the supernova rate in the underlying disc a more cussion of average properties of the mid-plane volume densities, 
or less established disc-halo interface is expected. We find that the gas distribution within the Milky Way H i disk 
Most prominent, comparing the gaseous content of galax- approximately exponential for a large part of the disk. The 
ies with the stellar distribution, is the large radial extent of the properties at > 35 kpc differ significandy; this is stud- 
gaseous disk. The outer part is only slightly affected by the stel- ^^d in Sect. 5. At the outskirts we find that the H i distribution 
lar disk but traces environmental influences such as influences "^eds to be described as a highly turbulent phase. In Sect. 6 we 
from the intergalactic medium. We explore global parameters consider the H i distribution at the polar caps. We require hydro- 
for layers above the disk and gas at the outskirts of the Milky static equilibrium and derive a velocity dispersion of <r = 74 
Way. For such an analysis we need to determine the properties km s 'for the halo gas phase. Such a phase explains not only the 
of the dominant part of the H i distribution first, i.e. the H i disk « ■ distribution at the outskirts but also the average local emis- 
which co-rotates with the stellar disk. Knowing the properties of ^lon high above the disk. About 10% of the H i is in this phase 

and we demonstrate that the Milky Way would look very similar 
to NGC 891 or NGC 2403 if it were observed edge-on at large 
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2. Database and derived density distribution 

Our analysis is b ased on the L e iden/A rgentine/Bonn (LAB) 
H I line survey (iKalberla et alJ l2005h . This survey com- 
bines the southern sky su rvey of the Instit uto Argentino 
de Radioastronomia (lAR) (iBaiaia et alJ l2005h with an im- 
proved version of the Leiden/Dwingeloo Survey (LDS) 
jHartmann & Burtonll 19971) . Currently, this is the most sensitive 
Milky Way H i line survey with the most extensive coverage both 
spatially and kinematically. The line profiles have been corrected 
for spurious sidelobe emission; we expect therefore that profile 
wings, most critical for our investigations, are only mildly con- 
taminated by instrumental problems. 

2. 1 . The density distribution n{R, z, <p) 

All the individual steps necessary to convert the observed bright- 
ness temperature distribution TbH, b, v) as a function of galactic 
longitude, latitude and LSR velo c ity to densities n(R, z, <p) have 
been described bv IKalberla et al.l (l2007h . We use cylindrical co- 
ordinates R, z, (p to describe the Milky Way disk; = is in di- 
rection / = 0. We use the lAU recommendations for the galactic 
constants; Rq = 8.5 kpc and v© = 220 km s"' . 

Figure [1] shows several global properties of the Galactic 
H I density distribution. We have chosen a common logarithmic 
transfer function to cover the high dynamic range of the density 
fluctuations within the Galactic disk. The average extension of 
the flaring disk layer is indicated. The thin white lines indicate 
the thi ckness of the disk at a 3cr level of the second moment fit- 
ted bv 'Kalbe rla et all ( |2007| ) . The warp is obvious for R > \5 
kpc as a sine-wave pattern. With increasing distance, the south- 
ern part of the disk, -180° 5 (p ^ 0°, returns to z ~ 0° while 
the northern part, Q° % ^ % 180°, bends up strongly and be- 
comes increasingly disrupted. Within the galactic plane several 
shells are detectable; some of them have very low central densi- 
ties. All kind of spurs, arcs, filaments, and chimneys are visible, 
indicating violent interactions for most parts of the disk, most 
pronounced for R ^ 20 kpc. 

We do not discuss details visible in Fig. [T]but a few com- 
ments concerning the TbH, b, v) to n(R, z, (p) conversion are nec- 
essary. The shape of the MiUcy Way rotation curve, but also the 
location of the Sun within the Galactic plane, causes complex 
geometrical transformations. A white line is plotted in Fig.[T]at 
latitudes \b\ - 30° +0.5° for visualization. Some ambiguities are 
discussed next. 

Local gas at low velocities may cause spurious features in 
the density distribution, most serious for |vlsrI < 10 kms"^ 
These regions have been disregarded, indicated by the black 
color in Fig.[T] Figure |2]displays the velocity field at mid-plane. 
The central green region, enclosed by isophotes, corresponds to 
I^lsrI < 10 kms"'. There are some more regions with unreli- 
able kinematic distances. In the direction of the Galactic cen- 
ter, we discarded data originating from \R\ < 3.5 kpc. This re- 
gion, approximately the extent of the bar, is affected by ambi- 
guities due to large veloc ity dispersions of the emission lines 
dWeiner & Sellwoodll 19991 Fig. 8). Low dv/dR values, predom- 
inantly at 1^1 ^ 30°, 4> <: 150°, and (p < -150° introduce ad- 
ditional distance uncertainties. We strictly exclude these regions 
from further analysis and use for the northern part of the Galactic 
plane flie range 30° < (p < 150° and 210° < < 330° for flie 
southern part. 

A fraction of the observed Hi gas distribution is presum- 
ably not associated with the disk. Intermediate and high velocity 
clouds (IVCs and HVCs, respectively) are most probably located 



above the disk dvan Woerden et al.ll2004l) . One might think that 
such features should be excluded from the database prior to the 
determination of n{R, z, (p) for the H i disk. We emphasize that no 
such filtering was applied to the observations. 

The scale height of the H i disk and the mid-plane density 
was determined by adaptive filtering of the volume density dis- 
tribution. The basic underlying assumption is that the disk can be 
characterized as a ho mogeneous ga s layer in circular rotations, 
for details we refer to lKalberla et d] (|2007). The thin white lines 
in Fig. [1] characterize the extension of the gas in the Galactic 
plane. We display the 3cr level of the second moment fitted by 
Kal berla et al.l (l2007l) to the flaring gas layer The 3cr limit is 
somewhat arbitrary, a different measure might be chosen to dis- 
tinguish between disk and extra-planar gas. Whatever criteria 
might be applied. Fig. [T] shows that most of the Hi gas settles 
down in a thin disk phase while there remains an extra-planar 
contribution, either diffuse or filamentary, but in most cases un- 
ambiguously associated with the disk. This extra-planar part was 
excluded when calculating disk scale heights and the mid-plane 
densities, but was left in the database for the rest of the analy- 
sis. There is direct evidence that on average 5% of the Hi gas 
is located outside the disk (Fig.[T2li. As discussed later, this is a 
lower limit only. 

z: O.OOOQOOe+OQ 
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Fig. 2. Milky Way velo city field at mid- plane for the mass distri- 
bution derived by Kal berla et al.l (|2007|) . The isophotes indicate 
velocities -190 < v^r < 190 kms~'in steps of 10 kms"'. The 
thick lines are at \visr\ - 50, 100 and 150 kms"'. 

Inside the solar circle the translation between distance and 
velocity is ambiguous because there are two positions along the 
line of sight that have the same velocity. These positions, how- 
ever, may have different z distances, so their densities may differ 
accordingly. We modeled the expected flaring of the H i distri- 
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Fig. 3. Derived mean surface densities of the Hi gas perpen- 
dicular to the disk (red). The dotted blue line gives Nh{R) in 
the northern, the dotted green line in the southern part of the 
Milky Way. The dashed black line repro duces the distribution 
published by 'Binnev & Merrifield| 119981 BM98) in Fig. 9.19. 
Wolfire et all (2003) (th ick black line) used partial data from 
Wouterloot etan(ll990l W90). 



bution and weighted the derived quantities correspondingly . Our 
method is comparable to that of Nakanishi & Sofud (|2003|) . ex- 
cept that they use a sech^ approach ( Spitzer 1942) in place of our 
more detailed model. Our model matches well wit h the observed 
flaring ; for a verification of the model we refer to lKalberla et alj 
(120071 Fig. 14). 

For a description of the global properties of the Galactic 
H I disk we discuss in the next sections the first three moments 
of the n(R, z, (p) distribution; the surface density, E(/?, 0), the 
first moment zq{R, <p), representing the mid-plane, and the sec- 
ond moment cr{R, (p), representing the scale height of the gas. 
In addition we derived characteristic mid-plane volume den- 
sities «()(/?, 0) for positions zoiR, </>), averaging over zo(^, <^) ± 
0.6 o-(R, 4>). 



3. The average radial surface density distribution 

Figure [3] displays average surface densities determined for the 
total disk, but also for azimuth 30° < (p < 150° (north) and 
210° < < 360° (south) separately. There are some deviations 
between north and south, we therefore derived the detailed sur- 
face density distribution across the disk in individual sectors of 

= 20°. Figure ID displays the results separately for both parts 
of the disk. Strong fluctuations, typically within a factor of two, 
with a correlation length of a few kpc are comm on. Some of 
the modulati ons are ca used by a spiral structure dLevine et alj 
l2006 a.b; Ka lberla et al.i.2007.) . 

Fig. |4] confirms the general systematic difi'erences between 
north and south. Most striking, however, are two regions with 
exceptional low fluctuations; in the north for \5 % R ^ 2Q kpc 
and in the south for 25 ^ 7? ^ 30 k pc. A comparison w ith the 
mass model and the flaring data from lKalberla et al] (12007 , Figs. 
18 & 19) indicates that these regions are located at the outer edge 
of the dark matter ring structures that were needed to explain the 
observed gaseous flaring. These structures have different loca- 
tions in the northern (/? ~ 13 kpc) and southern part (R ~ 18.5 
kpc) of the Galactic disk but extend in azimuth over at least 120°. 
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Fig. 4. Systematic variations of the average H i surface density 
distribution determined in 20° wide sectors for 30° < 4> < 150° 
(north) and 210° < (p < 330° (south). For comparison with Fig. 
[3]the fit (black) and total averages for north and south (blue and 
green respectively) are given. 

mean volume densities 
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Fig. 5. Derived average mid-plane volume density of the H i gas 
(red). The dotted blue line gives n{R) in the north (30° < (p < 
150°) , the dotted green line in the southern part of the Milky 
Way (210° < (p < 330°). The black line represents an exponen- 
tial fit with a radial scale length of 3. 15 kpc. 
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Fig. 6. Systematic variations of the average H i volume density 
distribution determined in 20° wide sectors for 30° < (p < 150° 
(north) and 210° < ^ < 330° (south). For comparison with Fig. 
|5]the fit (black) and total averages for north and south (blue and 
green respectively) are given. 



that the scale-height of the H i layer is constant. Some unspec- 
ified data were "arbitrary scaled by a factor of 2" by Dame to 
match othe r results. The case of a rising rotation curve was then 
adopted bv lBinnev & Merrifieldl d 19981 Fig. 9.19) as character- 
istic for the Milky Way. 

A slightly different surf ace density d i stribut ion, also plot- 
ted in Figl3] was used by IWolfire et al.] (|2003|) : we refer to 
their detailed discussio n. These au thors match 2=5 Mqpc"^ 
according to Dickev & LockmanI ([1990) with surface densi- 
ties derived by Wouterloot et al. ( Il990l) . Their curve rises for 
R ^ \Q kpc. At R ~ 12.5 kpc it matches our average sur- 
face density distri bution and at larg er distances our results agree 
well with iWouterloot et alj (Il990l) : their data are also plot- 
ted in Fig. [3] In gen eral we have fo r R ^ 10 kpc a goo d 
agreement wit h Hende rson et alj ( 1982h : Kulkarni et alj ( 1982 ); 
iDiplas & Savaga (il99ll) : IVoske'sldl999l) : lLevine et all (l200'6allbl) . 



Attempting to explain the discrepancies for ^ 10 kpc, the 
determination of H i surface densities close to RQmay be sub- 
ject to serious uncertainties, this was discussed in Sect. 2. In 
addition we found large biases caused by the assumption that 
the H I gas at |z| > kpc is stri ctly co-rotating with gas at the 
mid-plane (Kalbe rla et al.l l2007l see Figs. 3 & 5 for the flar- 
ing). Such biases may depend on the methods used. For a co- 
rotation we derive a 10% lower surface density for R % 12.5 
kpc, at the same time the exponential scale length increases from 
Rs = 3.75 kpc to Rs = 4 kpc. The s lope is consistent with 
previou s results (Henderson et al. 1982; Diplas & Savageiri991l : 
IVoskes|[r999: Levine et al..2006a.b) but the average surface den- 
sity for R ^ 12.5 kpc is still discrepa nt. We tested t he cas e of a 
rising rotation curve as proposed by I de Boer et alj (l2005l ) with 
a peak rotational velocity of 250 kms"^but found no signifi- 
cant changes of the local surface density. This most probably 
is caused by the fact that locally this rotation curve is flat. We 
question whether a rotation curve with a local rate dv/dR = 2 
km s^'kpc is realistic for the Milky Way. 



Despite all the scatter, the global radial surface density dis- 
tribution appears to follow an exponential law over a broad ra- 
dial range. We approximate the mean surface density distribu- 
tion by Z(/?) ~ 50 ■ e-(«-%)/f'^ wifli sq = 30 Mgpc^^ and 
Rs = 3.75 kpc. The exponential disk approximation is rea- 
sonable for 12.5 R ^ 30 kpc. For R ^ 12.5 kpc the sur- 
face densities saturate at approximately E,„„f,- ~ 10 M0pc"^. 
This result is in cl ear conflict with the conte mporary knowl- 
edge from reviews dDickev & Lockmanlll990l) and text books; 
we rep roduce in Fig. [3] the distribution of Binnev & Merrifieldl 
(119981 Fig. 9.19). Dickey and Lockman give 2=5 Mqpc"^, 
the average surface densities for R ^ 18 kpc according to 
iBinnev & Merrifieldl [[9m is even less. F or 7? ^ 18 kpc our 
data agree with lBinnev & Merrifieldl d 19981 Fig. 9.19); a general 
scale error cannot explain the problem. 

Discrepancies in Fig. [3] need a detailed discussi on, fortu- 
nately this is simplified greatly b y a recent review (LockmanI 
I2002al Sect. 4.3). Figure 9.19 of iBinnev & Merrifieldl d 19981) 
dates b ack to lDamd d 19931) who used data provided by LockmanI 
dl988 ^. Lockman studied two cases, a surface density distribu- 
tion for a flat rotation curve and alternatively the case of a curve 
that is flat for R <Rq, rising at a rate dv/dR - 2 km s^'kpc up 
to 2 Rq and staying flat afterwards. In both cases it was assumed 
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Fig. 7. Derived mean flaring of the H i gas (red). The dotted blue 
line holds for the northern, the dotted green line for the southern 
part of the Milky Way. For comparison an exponential fit for 
R < 35 kpc is plotted (long dashed) , also the flaring according to 



the model bv lKalberla et al.ld2007 
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4. Average radial mid-plane volume densities 

Figure |5] displays the average mid-plane volume density distri- 
bution n{R,ZQ) for the total disk, for 30° < 4> < 150° (north), 
and 210° < (p < 330° (south). Figure |6] shows how the volume 
densities depend on azimuth; these plots reproduce for R t 12.5 
kpc the same trends as seen before in the surface densities (Fig. 

For 1 % R % 35 kpc the average mid-plane volume density 
follows approximately an exponential distribution, n{R, zq) ~ no • 
^-(r-Rq)/r„ ^jjjj ^ Q 9 ^jjjj ^ 2.15 kpc. This line is 
drawn in Figs. |5]and|6] 

The surface density is essentially a product of mid-plane vol- 
ume density and the corresponding flaring. The fact that surface 
density as well as mid-plane volume density can be approxi- 
mated by an exponential function implies that the flaring also 
needs to be described well by an exponential function. We con- 
firm this as the fit /ir = /iq e'*"*©*/'^'' kpc with Iiq = 0.15kpc and 
/?o = 9.8 kpc for the half width at half maximum provides a good 
approximation to the observed flaring dKalberla et al.ll2007i, Fig. 
6) for 5 ^ R ^ 35 kpc. For R 35 kpc, data and fit diverge 
rapidly. F igure [7] re-plots for comparison the flaring derived in 
iKalberla et al.i (,2007,) . 

The finding that surface density, mid-plane density, and flar- 
ing can be represented by exponential functions does not reflect 
a general property of galactic disks but is closely related to the 
rotation curve used, hence it depends on the Galactic mass dis- 
tribution. In particular, the divergence of the exponential fit from 
the observed flaring for R ^ 35 kpc is caused by the shape of 
the gravitational potential which is oblate for small radii but 
changes to prolate in the outer part. The transition region, the 
location where we have a potential similar to the spheri cal case, 
is at R ~ 35 kpc (for a more detailed discussion see Kalberlal 
I2003L Sects. 5.2.2 & 5.2.3). This radius defines the "edge" of the 
disk. 



5. The Galactic outskirts 

Common to the average surface density and the average mid- 
plane density is that they both do not indicate a boundary of the 
Milky Way disk with a truncation of the observable H i distribu- 
tion. Beyond a radius of ~ 35 kpc the gas appears to have a more 
shallow distribution. We study this behavior in more detail. 
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Fig. 8. Average mid-plane volume densities for 20° wide radial 
sectors in the Milky Way. For comparison the fitted distribution 
of the HVD gas is given (red solid lines). 
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Fig. 9. Average mid-plane surface densities for 20° wide radial 
sectors in the Milky Way. For comparison the fitted distribution 
of the HVD gas is given (red solid lines). 



5. 1 . The warp, scalloping 

The Galactic warp was shown bv lKalberla et al.l (12007', Figs. 16 
& 17) to be well defined for R ^ 40 kpc. Three basic modes ex- 
ist with slowly variable phases. The outermost part of the warp 
for R ^ 30 kpc contains higher Fourier modes, thus the disk is 
scalloping. We determine these modes by subtracting the basic 
warp modes to 2 from the derived mid-plane positions. For the 
residuals we use a least square fitting algorithm to search for os- 
cillations with higher frequencies. We determine simultaneously 
amplitude, frequency, and phase. For various initial estimates in 
frequencies, corresponding to modes 3 to 20, we demand that 
the fit stays stable in frequency. As a result, we find oscillations 
with scalloping modes 10 to 20. The significance of these higher 
warp modes is low for R ^ 30 kpc. The associated power, de- 
fined by the squares of the amplitudes, is below 0.6% of the total 
warp power in this range. For R ^ 30 kpc we find a significant 
contribution; the preferred modes are 10 and 13. The relative 
strength of the scalloping increases linearly up to ~ 2.5% of the 



total warp power at = 40 kpc. This is shown in Fig. [TT] (ma- 
genta). The warp can be traced further with the scalloping power 
increasing to larger distances but the significance of the warp is 
highly questionable there. 

Our analysis is complementary to lLevine et al.l (l2006ah who 
used Morlet wavelet transforms for an analysis of the warp out 
to 7? = 30 kpc. They concluded that the m=10 and 15 scalloping 
modes are well above the noise, but localized; this suggests that 
the scalloping in this range does not pervade the whole disk, but 
only local regions. We applied a straightforward sine wave anal- 
ysis which is insensitive to local oscillations and found some co- 
herent grouping of higher mode oscillations, spreading out over 
several annuli, although insignificant for fits at individual radii. 
Only for R 30 kpc does scalloping become powerful enough 
to be significant in each case. This result supports Saha & Jo^ 
(|2006|) who concluded that scalloping at high modes has to be 
restricted to the outermost part of a galaxy. 
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5.2. Mid-plane and surface densities 

Figure |8] plots average mid-plane volume densities for 20° wide 
sectors as used before. Comparing the distributions with the ex- 
ponential fit from Fig.|5]we find a shallower decrease for R 35 
kpc. Average surface densities displayed in Fig.|9]show the same 
trend. Both plots suggest that the H i gas in the outskirts differs 
significantly from the Hi gas in the Galactic disk, which was 
fitted before by an exponential law. 

Tracing the shallow H i volume density distribution at large 
distances back to the observed brightness temperature distribu- 
tion we find that these belong to weak emission profile wings. 
The wings cannot be explained by emission from the cold 
or the warm neutral medium (CNM and WNM, respectively). 
Similar wings have been found in the Leiden/Dwingeloo Survey 
(Hartmann & Burton 1997) at high latitudes and Kalberl a et alj 
( Il998l Fig. 1) argued that the wings can be traced to low lati- 
tudes; this implies large distances. The LDS profile wings could 
be fitted well by a pervasive emission component with a high 
velocity dispersion (HVD) of ~ 60 km s ' . 

It is tempting to associate the shallow density distribu- 
tion, visible in Figs. [8] and |9l with such an HVD component. 
The proof, however, is not straightforward. The observed pro- 
file wings must originate from a diffuse distribution of unre- 
solved Hi clumps and it is a priori not obvious how these 
wings translate into the corresponding volume density distribu- 
tion. We therefore modeled the H i emission for an H VD com- 
ponen t with a radial scale length of Rhvd - 7.5 kpc (Kalberl^ 
l2003h . Synthetic emission lines were calculated for a dispersion 
of (Thvd = 60 kms^'and (Thvd - 74 km s"^ (differences will be 
discussed in Sect. 6.1). Next we derived the surface density and 
volume density distribution for the synthetic HVD distribution 
by applying the same procedures as used before for the analysis 
of the LAB observations. Figures [8] and |9] show the result. The 
HVD component is weak in comparison to the disk emission at 
low radii but dominates the H i distribution at the outskirts for 
R^35 kpc. 

The assumption of a pervasive extended HVD Hi compo- 
nent agrees well with the data (Fig. [8]and|9]l. The HVD model 
with (Thvd = 74 kms"' fits slightly better. As discussed before 
for the disk gas, there are systematic differences between north 
and south. Where is the HVD emission located in the observed 
Tb(1, b, v) distribution? 

In the Galactic outskirts the HVD gas becomes noticeable 
for R>35 kpc and the range 50° < < 130° and 250° ^ (f> ^ 
310° is useful for its analysis. The translation between Tb(1, b, v) 
and n{R,z,^), discussed in Sect. 2, is such that this gas must 
appear at high negative velocities for longitudes 40° $ / % 140° 
and at high positive velocities for 310° ^ / ^ 220°. Figure [TOl 
gives a few examples for the observed H i brightness temperature 
distribution, l-v slices at latitudes b - -10°, -5°, 5°, and b - 
10°. 

Extra-planar Hi gas has a low volume filling factor. 
Correspondingly, a patchy low level emission is expected in the 
region of interest. A second property is the high velocity dis- 
persion for this clumpy medium. This shows up as extensions 
(stripes) in the v direction. Figure[TO]shows a rather complicated 
pattern. A part of the emission may originate from large dis- 
tances but the warp makes it difficult to determine whether in- 
dividual filaments belong to the disk (having a low velocity dis- 
persion) or not. For a clearer picture we calculated the expected 
average disk emission and plotted the positions where this emis- 
sion should reach a level of 100 mK (red isophotes). Then we 



include the average emission according to the HVD model in 
Fig.|9](white isophotes). 

In the inner Galaxy the HVD emission is clear in Fig.fTOlwith 
high velocities against the disk emission at longitudes I ~ 25° 
and / ~ 335°. The region at I ~ 335° is currently studied by 
Alyson Ford with the Park es telescope. The regio n I ~ 25° 
was investigated in detail bv Lockmanl( l2002bl 12004). Numerous 
small clumps can be traced with the GBT up to distances of sev- 
eral hundred pc above the mid-plane. These clumps are usually 
unresolved in the LAB survey but contribute to the observed sur- 
face densities. 

Comparing the HVD model emission with observations, we 
find for some regions significant deviations from the expected 
distribution. This reflects the fact that HVCs and IVCs are 
grouped in complexes, while there are also regions with little 
extra-planar H i gas. 

5.3. Asymmetric gas distribution 

Figures [3] to |9] show a pronounced asymmetry of the gas distri- 
bution for ^ 20 kpc. For a more quantitative comparison we 
calculate the south/north ratios (S/N) for mid-plane volume den- 
sities, surface densities and for the scale heights. Fig.fTTIdemon- 
strates that the asymmetry is strongest for mid-plane volume 
densities. 
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Fig. 11. Asymmetry of mid-plane volume densities, surface den- 
sities, and scale heights characterized by their south/north ratios 
S/N. In addition the warp scalloping power in % is given. 

We argued in Sect. 4 that the H i gas in the Milky Way be- 
haves on average as expected for an exponential gaseous disk; 
the surface density is correlated with the product of the mid- 
plane density and the scale height. This balance, however, does 
not apply for the asymmetries visible in Fig.[TT] there is more 
gas in the south than in the north. The south/north asymmetry is 
even more pronounced if we compare the extra-planar gas frac- 
tion for R 25 kpc (Fig. [T2l l. The northern part of the extra- 
planar gas appears abruptly torn away to the south for R ^ 25 
kpc, in the direction of the Magellanic System. 

6. Local gas high above the disk 

The distribution of planar and extra-planar gas is locally aff'ected 
by dynamical processes. Here we consider a steady state that is 
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Fig. 12. Fraction of extra-planar gas observed outside the disk 
for -15 <z< 20kpc. 



dictated by the necessity of a global balance between gas pres- 
sure and gravitational forces. We discuss first how a pressure 
balance can be achieved by a multi-phase medium, next we use 
H I line data to verify the derived model. We also compare the 
extra-planar H i gas in the Milky Way with the anomalous gas 
phase observed in a few galaxies. 

6.1. Gravity and vertical gas pressure 

For a multi-phase medium, forces acting on the gas due to the 
pressure gradient p"^d p/dz need to balance the gravity (e.g. 
iBoulares & C"oxlll990l) and we may use this property to fit the 
parameters of the gas distribution . We use the mass distribution 
provided bv Kalb erla et al.l (12007 ) to calculate k. at R - 8.5 kpc 
up to z = 10 kpc (Fig. [T3T l. The volume density of the extra- 
planar gas phase is three orders of magnitude below the density 
of the other material in the disk. Correspondingly, k, is unaf- 
fected by the halo gas (at the peak, z ~ 5 kpc, the influence is 
^ 0.3%); the extra-planar gas is a genuine tracer of the Galactic 
potential. 

We fit p^^dp/dz for a multi-component gas to match k, 
and obtain agreement with a mean error of 0.6%. In Fig. [13] 
p^^dp/dz is over-plotted on k.; the gas pressure is corrected for a 
helium abundance of 9% of the hydrogen by number Essentially, 
both curves are identical, with differences hardly visible. Table 1 
lists the different phases, densities n and velocity dispersions 
(Tfi, used to balance the gas pressure gradient against the grav- 
itational field. For comparison we give nobs and (Tohs, derived 
previ ously by fittin g the gas distribution in the Solar neighbor- 
hood (lKalberlal200 3, Table 1). These parameters agree well. The 
most significant deviation is found for the velocity dispersion 
of the HVD component. This is a sensitive parameter in fitting 
the pressure gradient, but is less well determined when fitting 
observed profile wings. In Sect. 5.2 we used both solutions for 
comparison. 

The shape of k ^ for |z| $ I.l kpc is domin ated by the stel- 
lar distribution, see Kuiiken & Gilmord (Il989l) . We plot their k, 
relation for reference. At large z distances we find differences; 
the dark matter distribution dominates the shape of the k, func- 
tion and the gravitational forces become model dependent. The 
model assumptions by Kuijken & Gilmore (1989) lead to ever- 
rising forces which are physicall y implausible . The forces have 
to decrease at large z distances (lRohlfsllI977h . Such a solution 



was provided by iDehnen & BinnevI ( Il998l Galactic Model 2) 
and we plot their k, for reference. This model is currently fa- 
vored for ISM research (Cox 2005, Sect. 2.2) and the downward 
bending of k, is interpreted as a shaping provided by the non- 
pl anar geometry i n the Dehnen & Binney model. In the model 
of lKalberlaetal] (120071) the decrease of the forces at large z dis- 
tances is caused by the mass associated with the disk. The three 
k^ curves are discrepant for fe, <: 1.1 kpc. Stars cannot be used to 
decide between the models but we may use the gas distribution 
for this purpose. 




Kuijken and Gilmore (1989) 

dark matter disk with ring 

multi component gas 

no gaseous halo 

Dehnen and Binney (1998) 
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Fig. 13. Gravitational acceleration k. at R0 = 8.5 kpc, (solid 
red line) and gas pressure balance p ^^dpjdz (dotted green). For 
comparison we pl ot k^ according to iKuijken & Gilmord (Il989h 
(dashed, blue) and lDehnen & Binneyn[l998 l Galactic Model 2). 
The lower dotted line gives p^^dp/dz for a multi-phase gas with- 
out halo components. 



Attempting to fit the gas phase to match k- from 
IDehnen & BinnevI (1 19981 Galactic Model 2), we need only 80% 
of the observed density for the hot halo phase, but 5.5 times the 
observed density for the neutral halo phase. This, clearly, does 
not match the observations. We conclude that the observed extra- 
planar gas phase can be used to constrain dark matter models. 
Our approach is consistent with the extra-planar gas as observed 
in the local vicinity. However, some caution is needed. It is not 
expected that parameters derived at the position of the Sun agree 
exactly with model parameters for R ~ 8.5 kpc, since results 
are dependent on local influences, predominantly caused by the 
local hot bubble. A hydrostatic model can predict only global 
properties. 

Current ISM models (e.g. 'Coxl l2005l Sect. 2.1) do not in- 
clude a hot or neutral halo phase. For comparison, we plot 
p^^dp/dz for the fitted parameters from Table 1, but without a 
hot or neutral halo phase. This makes a huge difference. For a 
global pressure balance at large z distances a halo gas phase is 
unambiguously needed. Preliminary evidence for such a phase 
from in terstellar a bsorption li nes w as given by Miinch & Ziri^ 
(Il96lh . iLockman & GehmanI (1199 lb have shown that more than 
two Hi components are needed to explain the average local 
H I emission. In the next section we discuss the observational ev- 
idence for a multiphase Hi component in more detail. We show 
H I emission resulting from a neutral halo gas phase and demon- 
strate how the observations compare with the model given in 
Table 1. 
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Table 1. Local multi-component gas parameters used for Fig. 
[T3] in comparison to parameters derived from observations 
(Kalberla 2003) . 
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6.2. The polar caps in H i 
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Fig. 14. Average H i emission to logarithmic scale for latitudes 
b > 20°. The data are from the Bell Labs telescope (thick black 
line) and from the LAB survey (red). The step function shows the 
weights for each channel in dB. The thin red dotted line repre- 
sents broad Gaussian components indicating baseline problems. 
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Fig. 15. Average H i emission in logarithmic scale for latitudes 
\b\ > 20° from the LAB survey. The red line is for the north 
polar cap, the dashed green line for the southern part. The thin 
dotted black line represents the model according to Table 1 . 



Figure fT4l displays the average Hi emission measured at the 
north polar cap for b > 20°. We use logiTg) to compress the huge 
dynamic range of the^Hi emission. We compare data from the 
Bell labs survey (Stark et al. 1992, Fig. 10) with data from the 
LAB survey tKalberla et al. 2005). Both agree well, except for 
the wings at high velocities. Observational shortcomings cause 
the discrepancy. The bandwidth of the Bell labs filter bank was 
too small to cover the necessary velocity range completely; it 
needed to be shifted in velocity to obtain the main H i emission. 
The outer two channels of the filter bank were defined to be zero 
for all spectra. The weights (in dB, dashed line) show how this 
affects the mean Bell labs emission spectrum. As the weights 
decrease, the average Bell labs profile becomes incomplete and 
biased. 

The LAB survey may be influenced by systematic errors, 
too. The lower red dotted curve in Fig. [14] shows the mean of 
all broad LAB Gaussians, most probably representing baseline 
defects in the LAB survey (jiaiid & K alberla,2006) . This is two 
orders of magnitude below the mean emission; the weak profile 
wings are significant. 

Figure[T5]compares the LAB data for the northern and south- 
ern polar caps. We also plot the model distribution according 
to our best fit mass model, using the fitted parameters for the 
Hi components from Table 1. Note that the fit was to match 
p^^dp/dz with k-, but no fit was applied to the plotted profile 
in Fig. [15] and to the wings in Figs. [8] & [9] All profiles in Fig. 
[T5]show extended wings, but with some systematic differences. 
The enhanced emission for v;,,- ^ 100 km s ' at the southern po- 
lar cap is due to the Magellanic Stream. We find enhancements 
at visr ~ -60 kms"', caused by distinct IVCs and at visr ~ -120 
kms"', caused by HVCs. 

The broad wing underlying the CNM and WNM emission 
in Fig. [15] is caused by extra-planar gas, the HVD component. 
Numerous clumps, filaments, and spurs, located above the CNM 
and WNM disk, contribute to a signal which can be well repre- 
sented by a weak component with a large velocity dispersion. We 
conclude that the local observed extra-planar H i gas distribu- 
tion, forced to be on average in pressure balance k~ ~ p^^dp/dz 
(Fig. [TJt , is consistent with the gas component that was pro- 
posed in Sect. 5.2 to explain the gas distribution at the outskirts 
of the Galaxy. Systematic differences are caused by a warped 
and asymmetric gas distribution. Gas in the southern part of the 
outer Milky Way appears to be enhanced. This is consistent with 
the enhanced positive velocity wing in Fig.[T5]for b < 20°. The 
negative profile wing is enhanced for b > 20° due to the warp 
that bends up to the north. 

6.3. Extra-planar gas fraction 

During our analysis we distinguished disk emission and extra- 
planar gas. For a safe determination of the first moments of 
the stationary disk gas we masked and excluded outliers of the 
H I volume density distribution which are most probably not as- 
sociated with the disk. This part of the n(R, Z, (p) distribution has 
radial velocities that are incompatible with a circular rotation of 
the disk gas around the Galactic center Individual features show 
up as sheets and filanients in the n{R, Z, (p) distribution (see Fig. 
[Tlor iKalberlaetalJ dlOOTl Sects. 2.2 & 3 ). The coiTesponding 
clouds, visible in the Tb{1, b, v) distribution (Fig.[TOb. are usually 
classified as HVCs and IVCs. In general it is assumed that IVCs 
are located in the lower Galactic halo while HVCs are at larger 
distances (Albert & Danly 2004). 

For a determination of the extra-planar gas fraction we use a 
range where the properties of the Galactic disk are well defined. 
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10 ^ /? ^ 30 kpc. On average we find there direct evidence that 
about 5% of the Hi is not associated with the disk (Fig. [T2l i. 
This estimate is independent of the mass model used. Using pa- 
rameters from Table 1, integrating for R ^ 50 kpc, we derive 
that the HVD component amounts in total to 10% of the gas. 
This discrepancy is explained by the fact that most of the HVD 
gas phase in the Milky Way cannot be distinguished from the 
CNM/WNM gas phase. Only gas at the wings of the HVD dis- 
tribution is counted and we estimate that about 50% of this gas 
remains undetected by our filter and revise the extra-planar gas 
fraction to be 10%, consistent with the HVD model. 

6.4. Anomalous gas in the Milky Way, an external viewpoint 

To compare the Milky Way with sensitive observations of other 
galaxies, we take the viewpoint of an external observer. The 
H I emission is synthesized accordingly. We use a flat disk with- 
out warp at an inclination of 80° and calculate the total emission 
along the line of sight. The contours in Fig.[T6]represent observ- 
able column densities, normalized and plotted in steps of 2 dB. 
We consider observations with a limited dynamic range of 30 
dB; the lowest contour is at -30 dB. 



total HI emision 




5 10 15 20 25 30 
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Fig. 16. Position velocity diagram along the major axis, observ- 
able for a galaxy like the Milky Way with an inclination of 80°. 
The thin contours represent normalized column densities in steps 
of 2 dB; the lowest contour is at -30 dB. The green dashed lines 
gives the rotation curve at z = kpc, the solid red line the mean 
rotation for the H i in the halo. 

The simulated position velocity diagram shows that the ve- 
locity at the peak of the line emission is identical with the ro- 
tation curve from the model. The average observable rotation 
velocity for the halo gas phase alone, determined from the peak 
emission of the halo gas phase lags behind. The total column 
density distribution is symmetric at large column densities, but 
as a consequence of the lagging halo some asymmetries appear 
at low levels. Projection effects cause extended regions with 
Hi emission at lower velocities. Such features at levels below 
-20 dB were first obs erved in NGC 891 dSwaters et al.]|l997h . 
the so-called beards ( Fraternali et al. I l2001h . Close to the center, 
the beards extend to forbidden velocities, mimicking an apparent 
counter-rotation. 

Beards and forbidden velocities in Fig. [16] have the same 
origin. The line of sight passes through various layers of extra- 
planar gas that rotate slower than the disk. The lag is approxi- 



mately constant for zlK {R) with z distances normalized to scale 
heights (l Kalberlal2003L Fig. 11). For small radii, the densest part 
of the thick extra-planar gas disk close to the center becomes 
most prominent. The lag is most pronounced there. 

Our model takes both the flaring and lagging halo into 
account. Both are necessary to explain the observations. The 
beards in Fig. [16] originate from the same Hi halo gas phase 
that causes the extended wings in Figs. [14] & [15] At the end of 
this subsection we point out that an external observer, with a 
telescope that is limited in dynamic range to 30 dB, would see 
the Galactic Hi disk out to ~ 30 kpc only. At low levels it 
is far easier to find e xtended emission perp endicular to the disk 
as shown recently bv lOosterloo et al.l (12007 ') in the case of NGC 
891. Reaching the outskirts at distances of /? ~ 50 kpc, is a de- 
manding task and needs a dynamic range of 40 dB or better The 
HVD component at such radii is essentially invisible. 



7. Summary and discussion 

We used the LAB survey to derive the Hi volume density dis- 
tribution in the Milky Way. The Tb{U b, v) to n{R, z, (f>) conver- 
sion depends on the rotation curve; we have chosen an almost 
flat curve with lag ging rotation for gas above the mid-plane 
dKalberla et al.l 1200 7). The corresponding mass model implies 
that some of the dark matter is associated with a thick exponen- 
tial disk. 

Our aim was to derive global properties for the H i gas distri- 
bution in the Milky Way. We calculated the first moments to de- 
termine surface densities, mid-plane positions, associated mid- 
plane densities, and scale-heights. Planar and extra-planar gas 
was distinguished by applying adaptive filters. We discuss mid- 
plane positions first. 

7.1. The warp 

The Galactic warp can be traced from R ~ 12.5 kpc out to R ~ 
45 kpc but the coherent structure is lost for R k, 40 kpc. The 
scalloping power, defined by the squares of the amplitudes from 
a sine wave fit, is below 0.6% in the inner part of the disk. For 
30 ^ ^ 45 kpc we observe a linear increase up to 3% of the 
total warp power. 

7.2. The exponential H i disk 

We derive average surface densities and, taking the warp into 
account, average volume densities at mid-plane. In both cases 
we find significant fluctuations and asymmetries between north 
and south which are correlated. Surface densities and mid-plane 
volume densities are affected by spiral structure, but differences 
betw een north and south ne ed to be explained by an asymmetric 
disk dLevine et al.l l2006a"b'; 'Kalb erla et alj|2007l) . 

Average surface densities, mid-plane volume densities, and 
scale heights can be approximated over a broad radial range by 
exponential relations. For the surface density we find a scale 
length of 3.75 kpc within 12.5 ^ R < 30 kpc. Mid-plane vol- 
ume densities are well approximated with a scale length of 3.15 
kpc for 1 ^ R ^ 35 kpc. Both relations imply that the scale 
heights also can be described well by an exponential relation, in 
this case valid for 5 ^ R ^ 35 kpc. 

For R ^ 12.5 kpc we find approximately constant surface 
densities (J.i„„er ~ 10 Mqpc"^). This is also the inner part of the 
H I disk which is unaffected by the warp. 
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7.3. The Galactic outskirts 

The exponential disk approximation for the disk gas breaks 
down in all respects at radial distances R ^ 35 kpc. In the outer 
part we find a shallower distribution which can be approximated 
by a gas phase with a radial scale length of 7.5 kpc and a high 
velocity dispersion of 74 km s"'. 

Ou r results are apparently in disagreement with iMalonevI 
( 1199 3*). who found a sharp truncation of the neutral hydrogen 
distribution for the outer disk Hi gas in NGC 3198 at a surface 
density level 10"^ ^ Nh ^ 10'^ cm"^. This Hmit is attributed 
to the extragalactic radiation field. To study the origin of the 
HVD component at large distances we examined the data. At 
a distance of 60 kpc the telescope beam covers a typical area of 
600 pc in diameter Any emission within this area is smoothed 
out completely. Inspecting the data cube visually, we find re- 
gions with peak densities of n{R) ~ 10^^ cnT^, surrounded by 
less dense envelopes. These unresolved individual condensations 
may be well above the predicted cut-off limit. 

Maloney's interpr etation was questioned recently by 
lOosterloo et alj (l2007b who found that a deeper mapping of 
NGC891 led predominantly to the detection of unknown low 
level emission in the vertical direction, while the radial limits 
remained. Such an asymmetry is unexpected for an isotropic ra- 
diation field. Our analysis shows the same trend. At the outskirts 
of the disk the HVD component becomes most extended in z di- 
rection. Mid plane volume densities are dominated by the HVD 
component for R ^ 40 kpc, while surface densities (integrated in 
z direction) show the same effect for Rk. 35 kpc (compare Figs, 
landig. 



7.4. Polar caps and extra-planar gas 

The Hi gas distribution, averaged over all latitudes \b\ > 20°, 
probing approximately the local volume within a radius of 1 kpc, 
has broad profile wings. These wings may be affected by stray 
radiation of the telescope, but for the Bell Labs and the LAB 
survey they are not explainable by instrumental problems. The 
same HVD component with a dispersion of 74 kms which 
was used to fit the gas distribution in the outskirts, is useful also 
to explain this local component. A multi -phase medium contain- 
ing this HVD component is found to balance gravity by its pres- 
sure gradient, ~ p ^^dp/dz for the mass model proposed by 



iKalberla etafl(l2007h . 



7.5. Clumps, fast and at large z distances 



Evidence for a population 
locity dispersion of cr 
[Radhakrishnan & Srinivasai 
( 1985); Lockman & Gehmai 



of fast gas clouds 
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with a 
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[Kalberla et al. ( 1998), analyzing the Leiden/Dwingeloo Survey, 
proposed a more extreme dispersion of cr ~ 60 kms ' which 
is modified here to cr ~ 74 kms~'. Most remarkable, how ever, 
is the early prediction by iPikelner & Shklovskvl (Il958l) who 
proposed a neutral halo gas phase with cr ~ 70 km s ' . 

The HVD component analyzed here as a diffuse emission 
component with a dispersion of cr ~ 74 kms"', is patchy. It is 
made up of a large number of filaments and small Hi clumps. 
These stand out from the normal disk emission only if they 
reach large z distances or if they have high forbidden veloci- 
tie s. Evidenc e for su c h clumps in the inner Gal axy was given 
by iLockmanl (l2002bl) . IStanimirovic et al.] (l2006l) found similar 
clumps in the anti-center direction. Clouds with high forbidden 



velocities in the inner Galaxy were detected bv lStil et aTl (l2006h 
in the VLA Galactic Plane Survey. These authors estimate that 
this population of fast clouds has a radial scale length of 2.8 to 8 
kpc, consistent with the scale length of 7.5 kpc derived by us. 

7.6. Extra-planar gas fraction 

For our analysis we used adaptive filters to distinguish the sta- 
tionary disk component from extra-planar gas. Over the main 
body of the disk (10 ^ ^ 30 kpc) we measure that on aver- 
age 5% of the gas is outside the disk (Fig. [T2]). This is a lower 
limit only. Assuming that the extra-planar gas can be character- 
ized as an HVD component implies that a significant fraction of 
this phase is hidden from direct observation since it cannot be 
distinguished from disk gas. For our analysis we estimate that 
only 50% of the LVD gas can be discriminated, resulting in an 
extra-planar gas fraction of 10%. This is consistent with a hydro- 
static model (Sect. 6.1), observations at the polar caps (Sect. 6.2) 
and also with the gas distribution at the Galactic outskirts (Sect. 
5.2). We conclude that all over the Galactic H i disk the "normal" 
phase, the CNM and WNM, is associated by an "anomalous" gas 
phase, described here as an HVD phase. For the inner Galaxy 
(/? $ 10 kpc) we observe an increase of the extra-planar gas 
fraction of a few percent. This might be explained by fountain 
activities but, as discussed in Sect. 2, estimates at such distances 
are highly uncertain. 

In this light, the general features of the extra-planar H i gas 
in NGC 891 and NGC 2403 appear to be comparable (within a 
factor of 2) to what we see in the Milky Way. This is confirmed 
by projecting the Hi gas distribution of the Milky Way to large 
distances. Simulating position-velocity diagrams, the HVD com- 
ponent causes wings and beards close to the center. The gaseous 
outskirts, however, would hardly be observable. The telescope 
sensitivity needs to be improved by a factor of 10 to observe this 
part, which is feasible with the SKA. 

7.7. HVCandlVCgas 

HVCs and IVCs are usually classified from observations accord- 
ing to their velocities. One can distinguish between coherent 
complexes and isolated compact objects. Their origin may be 
from fountain flow or accretion. Also, metalicities are important 
but here we discuss H i properties only. 

For our analysis of the global properties of the H i distribu- 
tion we have chosen different criteria. The primary goal was to 
derive the properties of the main body, the H i disk. This has a 
well defined dynamically cool gas distribution. There are two 
phases, CNM and WNM, both with velocity dispersions perpen- 
dicular to the disk that are low compared to the rotation velocity. 
The gas disk is flaring but the most prominent inner part of the 
H I disk is highly flattened. 

However, as proposed bv lSpitzeJ (11956 *) and much later es- 
tablished by X-ray observations of the predominantly hot ISM 
of elliptical galaxies, the hot Galactic ISM does not have to be 
so flat. The extra-planar H i is embedded in the hot component. 
This arises the question of whether there might be a quasi-static 
system near equilibrium as modeled by Spitzer or whether the 
Milky Way halo is strongly affected by steadily changing condi- 
tions, e.g. due to mass inflow or star bursts. 

The star-formation-rate (SFR) in the solar neighborhood is 
remarkably constant over the lifetime of the disk, infall of mate- 
rial a ppears also to b e approximately constant at one-half of the 
SFR (lTwarog|ll980l) . Ongoing accretion may be derived from 
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HVC complexes and also from the Magellanic Stream but it 
is unclear whether irregularities in accretion rate and outflows 
are strong enough to aff'ect the global shape of the Galactic halo 
significantly. Little is known about the impact of star bursts on 
the structure and dynamics of the Milky Way halo. Supershells 
and chimneys may reach z distances of several kpc but ap- 
pear not to be able to shape a significant fraction of the halo 
CMcClure- Grifliths et al. 2006|). 

As a halo model therefore we have chosen the simple case 
of a quasi-static system that can be approximated numerically 
by a hydrostatic distribution. We are interested in global prop- 
erties of the gas distribution, not in details that are more im- 
portant for individual components. A global model makes sense 
when averaging the H i emission over several kpc as we obtain a 
self-consistent solution. As discussed in Sects. 5 & 6, the same 
model parameters apply to both polar caps and to the sectors 
5Q° <4>% 130° and 250° ^ < 310° at40 > > 60 kpc. These 
are the only regions available for an unambiguous separation of 
planar and extra-planar gas in the Milky Way. 

Figure [TO] a l-v plot of measured brightness temperatures, 
displays numerous IVCs and HVCs. These features appear 
partly to be organized in complexes but are in general close to 
the disk. Projection effects cause a complicated picture, in par- 
ticular velocity crowding in some parts of the sky makes it dif- 
ficult to judge whether a particular feature belongs to the disk 
or not. Deviation velocities between planar and extra-planar gas 
components depend strongly on position. 

The density distribution n{R,z,(p) in Fig. [1] shows a clearer 
picture. Extra-planar gas is offset from the Galactic disk but most 
of it is clearly associated with the disk. Again we find large scale 
fluctuations but there is a general trend that the extra-planar gas 
distribution increases with decreasing mid-plane distance. This 
situation appears best to be described as a "thick disk" having 
twice the radial scale length of the CNM/WNM gas distribu- 
tion. We use the model of a turbulent HVD distribution. Local 
enhancements, or complexes, are considered as statistical fluc- 
tuations; they become less dominant if one averages over larger 
areas, this is a typical property of a turbulent medium. FigurefTSl 
showing average emission of the polar caps at latitudes \b\ > 20° 
also includes IVCs, HVCs, and the Magellanic stream. These 
local phenomena show up as minor perturbations. 

Comparing properties of the HVD component with global 
properties of the HVC gas we find some interesting agreements. 
HVCs have a multi-phase structure an d are most probably con - 
fined by an external hot halo phase (iKalberla & HaudI l2006h . 
HVCs have specific turbulent energy densities that are an or- 
der of magnitude higher than that of comparable clumps in the 
Galactic disk. The same conditions hold for the HVD compo- 
nent. Most surprising, however, is that the integrat ed flu x distri- 
bution for HVCs as function of velocity (Wa kkeJI 19911 Fig. 5) 
is a continuation of the gas distribution in the broad wings of the 
HVD component to high velocities (Fig.fTsTl. 

The HVD component is coupled to the disk but it is anoma- 
lous in the sense that it represents a highly turbulent clumpy 
phase, best described as an extended very thick disk. At the po- 
sition of the Sun we determine an exponential scale height of 3.9 
kpc. The HVD component can be best characterized as an extra- 
planar gas component; it expands the phase space occupied by 
the H I disk predominantly perpen dicular to the disk. IV Cs are 
located in the lower Galactic halo (lAlbert & Danlvll2004 . Their 
velocities as well as known z distances place them at the less 
distant part of the HVD phase space. Distances for HVCs are 
less well known but according to their velocities at least some of 
them may belong to the outer part of the HVD phase space. The 



HVC velocity limit is \visr\ - 465 kms roughly coincident 
with the escape velocity of the Milky Way. 

Individual extra-planar gas layers are extremely patchy and 
lag behind the underlying rotating disk. These layers are aff'ected 
by flaring and the lag depends on the distance zjh^iR), normal- 
ized to the disk scale height. The anomalous gas phase is rather 
extended. We could trace it up to /? ~ 60 kpc. 

The gravitational potential at the outskirts of the Milky Way 
is very shallow. The rotation of anomalous gas for R ^ 35 kpc 
is therefore essentially unaff'ected by flaring. The fact that we 
observe an extended HVD gas phase at 40 ^ ^ 60 kpc implies 
accretion by a "drizzle" of low mass fragments. This gas cannot 
originate from fountain flows. Fountain events far outside R25 
are highly improbable as the radial scale lengths for stars and 
HVD gas differ significantly. 

A continuous flow of material i n galaxies f ro m th e 
disk into the halo was studied by Collins et"an ( |2002|) ; 
iFraternali & BinnevI (|2006'); 'Hea ld et all (|2006a'). Their models 
can explain the vertical gas distribution but fail to reproduce the 
apparent counter rotation and vertical velocity gradients. They 
also produce a general outflow, contrary to observations. 

7.8. The lopsided H 1 distribution 

Our aim was to describe global features of the Hi distribution 
that may be representative of a steady state. Calculating radial 
distribution functions we found in all cases significant system- 
atic differences between the northern and southern part of the 
Milky Way. The system is asymmetrical. Midplane volume den- 
sities and surface densities are most discrepant at /? ~ 20 kpc. 
Considering the Hi flaring we found a systematic modulation 
that may be explainable by a misaligned gravitational field orig- 
inating from a misaligned dark matter disk (Fig. [TTT i. The flar- 
ing ratio peaks a t R ~ 25 kpc are in excellent agreement with 
IWeinbera (1 19981) who predicted that the interaction between 
Milky Way and Magellanic system causes a dark matter wake 
at this position. 

7.9. Rotation curve and mass distribution 

Our results depend on the rotation curve, hence on the mass 
model. We tested several different possibilities and summarize 
briefly how much our conclusions are affected by the mass 
model used (Kalberla et al. 2007). 

The global surface density distribution, mid plane d ensities 
and also flaring depend weakly on the rotation curve (Voskes] 
1999; Nakanishi & Sofue 2003; Levine et al. 2006a b) but our 
results are not altered significantly as long as a nearly flat rota- 
tion curve is used. The radial exponential scale length of 3.75 
kpc derived by us changes to 4 kpc if one uses a constant rota- 
tion with a circular velocity of 220 kms"'. Our conclusion that 
surface densities, mid plane densities and flaring can be approx- 
imated by exponential relations remains valid, although the cor- 
relation degrades somewhat if we use alternative rotation curves. 

Deviations from axisymmetry depend on non-circular terms 
in the rotation curve. We used the first order epicyclic streamline 
correction bv iLevine et al.l (l2006a . Eqns. 1&2) which is a great 
improvement in overcoming discontinuities at large radial dis- 
tances. This correction is empirical but the assumption that sur- 
face densities should not have azimuthal discontinuities remains 
valid for all mass models. The correction is consistent with an 
asymmetrical mass distribution and the results do not depend 
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significantly on the details of the mass model. The same holds 
for asymmetries in the H i distribution. 

The average Hi emission at the polar caps (Sect 6.2) has 
broad high velocity wings. This is known from observtions. 
Parameters for the underlying HVD component are completely 
independent of any mass model. The parameter fitting, assuming 
halo gas in global hydrostatic equilibrium with the gravitational 
field, depends strongly on the mass model. A dependence on the 
mass model exists also for HVD gas at the outskirts of the Milky 
Way. It is the shape of the outer rotation curve that determines 
the distance of this component. 

We argued that the thin Hi disk, defined by the CNM & 
WNM gas, is associated with a thick extended HVD gas disk. 
This model is self-consistent for the mass model used but it nec- 
essarily has to break down for significant changes of the mass 
model or the outer rotation curve. 
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Fig. 1. Derived volume density distribution at R - 10, 12.5, 15, 17.5, 20, 22.5, 25, 30 kpc (top to bottom). Locations affected by 
local emission with \visr\ < 10 kms"' are flagged dark. We display densities n < 0.5 cm"^. A logarithmic transfer function was 
chosen to emphasize low densities. The thin white lines indicate the average flaring at a 3cr level and enclose the warped disk. The 
thick white lines indicate latitudes \b\ - 30° + 0.5°. 
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Fig. 10. Position-velocity diagrams for the Milky Way Hi brightness termerature at latitudes -10°, -5°(left) and 10°, 5°(right). The 
white isophotes display HVD emission from the model (10% to 90% in steps of 10%), the red isophote marks an average disk 
emission of lOOmK. A logarithmic transfer function was used for Tb < 30K. 



